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Abstract 

Multi-fractal scaling in the transition to the dissipative regime for fully-developed com- 
pressible turbulence is considered. The multi-fractal power law scaling behavior breaks down 
for very small length scales thanks to viscous effects. However, the effect of compressibility 
is found to extend the single- scaling multi-fractal regime further into the dissipative range. 
In the ultimate compressibility limit, thanks to the shock waves which are the appropri- 
ate dissipative structures, the single-scaling regime is found to extend indeed all the way 
into the full viscous regime. This result appears to be consistent with the physical fact 
that vortices stretch stronger in a compressible fluid hence postponing viscous intervention. 
The consequent generation of enhanced velocity gradients in a compressible fluid appears to 
provide an underlying physical basis for the previous results indicating that fully-developed 
compressible turbulence is effectively more dissipative than its incompressible counterpart. 
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1 Introduction 



Compressibility effects on fully-developed turbulence (FDT) are of importance in modern 
technological flow processes like makes of supersonic projectiles, hypersonic re-entry vehicles 
and high-speed internal flows in gas-turbine engines. Astrophysical processes like star for- 
mation in self-gravitating dense interstellar gas clouds via Jeans' instability (Chandrasekhar 
[Ij and Shu [2j) are other cases in point. So, considerable effort has been directed on this 
problem to date (Moiseev et al. [3], Passot and Pouquet [5], Porter et al. [6], Erlebacher 
et al. [7|, Shivamoggi [8]- [15] and Kritsuk et al. [IS]- [IT]). 

The multi- fractal model was used by Shivamoggi [9]-[II], [H], [I5] to describe the nu- 
merically observed spatial intermittency in fully-developed compressible turbulence (Passot 
and Pouquet [1], [S] and Lee et al. [IB])- A multi- fractal object exhibits a global scaling 
structure that is described by a continuous spectrum of scaling exponents a and in a certain 
range / = [ttmm, OLmax\- Each a G / has the support set S{a) C of fractal dimension /(a) 
such that the velocity increment over a small distance I has the scaling behavior (Frisch and 
Parisi ^19]), 

bv{t) ~ r, femall. (1) 

The sets S{q) are nested so that S'(a') C S{a) for a' < a. The fractal dimension /(a) is 
obtained via a Legendre transformation of the scaling exponent of the pth-order structure 
function, 

SM^ <IH0\'> ^/..Hf— .^^^ (2) 

according to 

= inf [ap + 3 — /(a)], say for a = a*. (3) 

a 

where, 



p. (4) 



da 

The factor (3 — /(a)) is the fraction of the coarse-graining measure of the volume of the set 
S{a) over the total volume. 

When viscous dissipative effects arise, the multi-fractal power law described by equation 
(2) breaks down for very small £ and the cut-offs are determined by external parameters like 
the Reynolds number. In such cases, as in incompressible turbulence (Wu et al. [20], Frisch et 
al. [21]), the multi- fractal in question turns out to exhibit a pseudo-algebraic behavior with 
certain universal features, albeit on using a suitable rescaling via a multi-scaling method (Wu 
et al. [SO])- Similar situation arises also for a passive scalar diffusing in a random velocity 
fleld (Shivamoggi [22]). 



2 Multi-fractal Scaling in the Inertial Range: Behavior in the 
Ultimate Compressibility Limit 

In order to relate the fractal dimension f{a) to the generalized fractal dimension Dg of the 
kinetic energy dissipation fleld i, one considers a coarse-grained probability measure given 



2 



by the total kinetic energy dissipation ocurring in a box size 



(5) 



and then covers the support of the measure e with boxes of size £ and sums the moments 
over all boxes. Noting the asymptotic scaling behavior of these moments (Halsey et 
al. [23]) 

i=l 

or 



{q-l)D, 



(6a) 



dfi(a 



(6b) 



7 is the ratio of specific heats of the fluid and may be used as a compressibility parameter 
(1 < 7 < oo) - the incompressible fluid corresponds to the limit 7 ^ oo while the ultimate 
compressibility case given by the limit 7^1 corresponds to Burgers turbulence. 

The dominant terms in the integral in (6)b, in the limit of small £, may again be extracted 
using the method of steepest descent: 



37 



7 



a* + 2 



with, 



dfia*) 



q-f{a*) = {q-l)D, 

37- r 



da 



7 



(7) 



(8) 



Eliminating /(a) and putting q 

mi 

•7-r 



7-1 
37- 



Y ) we obtain from (3) and (7) (Shivamoggi [8] 



37 



^-3 



37-3 



p — 3 



(9) 



.37-1 

On noting the scaling behavior the density p (Shivamoggi [8]) (which follows on applying 
scale-invariance arguments directly to the Navier-Stokes equations for a compressible fluid 
in conjunction with the scale-invariance condition on the kinetic energy dissipation field), 



p ~ (5i;)T-i 

the scaling behavior of the kinetic energy is given by 

< p{6v)'^ > < 6v >^ 
(11), in conjunction with (2), leads to 



< p{6vf > ~ /(^) 



where, we have, from (9), 

'^(7^) 37 



27 



7-1 
37-1 



3 — D( 27 ' 



(10) 



'11^ 



(12) 



(13) 
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In the ultimate compressible limit 7^1, (12) yields 

< p{5vf > ~ £ (14) 
which leads to the Kadomtsev-Petviashvili [21] spectral law for compressible turbulence - 

E{k) ~ A;-2. (15) 
On the other hand, for a fractally homogeneous turbulence, we have 

Dg = Do,yq (16) 

so (9) becomes 

/ ^/ _ 1 \ 1 r /^'^o/ — ^'^N 1 

(3 -Do). (17) 



37 - 3 , „ 

p — 6 



(37-1)^ 3LV37-1, 

If we identify the dissipative structures in compressible turbulence with shock waves, we 
have Dq = 2, and (17) yields 

Cp = h^P- (18) 

(12), in conjunction with (18), leads to 

< p{5vf > ^ e (19) 

in agreement with (14). 

3 Multi-fractal Scaling in the Transition to the Dissipative Regime 

The dissipative effects materialize when typically the eddy turn-over time exceeds the viscous 
diffusion time, i.e., 

dv p 

or 

6v < ^ = u^ r. (21) 

(21) implies the presence of a viscous cut-off u so that the boxes with measures 6v < u ot 
indices a > a, where 

inu ^ ina , , 

would be empty and not contribute therefore to the integral for in (2). Note that different 
length scales (below the Kolmogorov microscale 77) have different viscous cut-offs. Thus, 

Spii,u)= j rf^(a)rP+='-^(") (23) 

which signifies the existence of a continuous range of viscous cut-offs in the multi-fractal 
model for the inertial range. As these cut-offs vary, a cross-over in the scaling of the structure 
function appears for ^ ~ 77 (Wu et al. [20]). 
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For compressible turbulence, noting that 



we have, 



S 2 

^— (24) 

ri 1]^ 



^ and ~ ( ^ ) ' (25) 



pv \p 
from which, 

(26) 

On the other hand, on noting the following scaling behavior (Shivamoggi [S]), 

p{l) ~ (27) 

and using (26), we have 

«^>,> or-^>^. (28) 



Using (28), we have from (22), 



7 + 1 

>r]: a> > a* (29) 

37 — 1 



where a* is defined in equation (3) and gives the leading order contribution to the scaling 
exponent of the structure function. So, in this interval, we have the usual single-scaling 
multi-fractal regime in the inertial range. On noting that 

7 + 1 1 I 

- + (30) 



87 — 1 3 87 

we observe, from (29), that the effect of compressibility is to extend the single-scaling multi- 
fractal regime further into the dissipative range, provided of course a^, < a. 
Further, in the ultimate compressibility limit 7^1, 

^ + ' - 1 (81) 



87-1 

so, according to (29), at the onset of viscous effects, we have 

dvr^i (32) 

which is, however, the scaling signifying a smooth velocity field appropriate for the full 
viscous regime. So, in the ultimate compressibility limit, thanks to the shock waves which 
are the appropriate dissipative structures, the single- scaling regime extends indeed all the 
way into the full viscous regime. 

On the other hand, when i falls well below r] so that a < a* < (^^^), the minimum of 
is attained for a = a. We now have a pseudo-algebraic behavior (as in the incompressible 
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case (Jensen et al. [25])) induced by dissipative effects - the power law in question has a 
scale- dependent exponent, albeit slowly varying via proportionality to ^n^. 
Thus, for small we have 

>p , if a > a,, 

S,{^,u)^{ (33) 
^sp+3-/(a) ^i{a<a^. 

It may be noted that the saddle-point evaluation method used in deriving (33) becomes 
exact in the limit a =^ otmin when one obtains the result corresponding to a homogeneous 
fractal. 

It is possible to bring out universal features even in a non-power-law situation, albeit on 
using a suitable rescaling associated with multi-scaling reflecting variations in the cut-off 
parameter u by introducing 

CnS fnf 1 

(33) then implies 

Fp{e) = I (35) 

[p + e[3-f{l)], if^>^. 

The first regime in (35) corresponds to the single-scaling regime where the graph of Fp vs. 9 
is a straight line with slope ^p. The second regime in (35) has the structure function data for 
different Reynolds numbers collapsing onto a single curve in the neighborhood of the viscous 
regime - multi-fractal universality for compressible turbulence. 



4 Discussion 



In fully- developed compressible turbulence, when the dissipative effects arise, the multi- 
fractal power-law scaling behavior breaks down for very small length scales and the cut-offs 
are determined by external parameters like the Reynolds number. However, if one uses 
a suitable rescaling, the multi-fractal in question is shown still to exhibit a multi-scaling 
behavior with certain universal features, as in the incompressible case, via the existence 
of a continuous range of viscous cut-offs. The effect of compressibility is found to extend 
the single- scaling multi-fractal regime further into the dissipative range - it is as though 
compressibility allows length scales to become smaller without forcing viscous dissipative 
effects to become operational. In the ultimate compressibility limit, thanks to the shock 
waves which are the appropriate dissipative structures, the single-scaling regime is found to 
extend indeed all the way into the full viscous regime. This result appears to be consistent 
with the fact that vortices stretch stronger in a compressible fluid hence postponing viscous 
intervention. The consequent generation of enhanced velocity gradients in a compressible 
fluid appears to provide an underlying physical basis for the results (via group-theoretical 
arguments applied to a Hopf-type functional equation formulation (Moiseev et al. |.3J), 
scahng arguments applied directly to the compressible Navier-Stokes equations (Shivamoggi 
[H], [S]), direct numerical simulations (Passot and Pouquet fSj, Erlebacher et al. and 
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Kritsuk et al. [16], [T7]) indicating that fully-developed compressible turbulence has steeper 
energy spectra and is effectively more dissipative than its incompressible counterpart. On 
the other hand, the enhanced vortex stretching is consistent with the results that - 

• the velocity-space singularities are stronger in fully-developed compressible turbulence; 

• the mean kinetic energy of a typical spectral mode in an equilibrium distribution is 
reduced by compressibility effects (Shivamoggi p!2]). 
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